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Abstract 

In this paper, we are concerned with optimal decay rates for higher order spatial deriva¬ 
tives of classical solutions to the full compressible MHD equations in three dimensional 
whole space. If the initial perturbation are small in i7^-norm and bounded in L'^{q G [ij f) )- 
norm, we apply the Fourier splitting method by Schonbek [Arch.Rational Mech. Anal. 88 
(1985)] to establish optimal decay rates for the second order spatial derivatives of solutions 
and the third order spatial derivatives of magnetic field in L^-norm. These results improve 
the work of Pu and Guo [Z. Angew. Math. Phys. 64 (2013) 519-538]. 
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1 Introduction 

In this paper, we are concerned with the compressible viscous and heat-conductive magneto¬ 
hydrodynamic (In short, MHD) equations in the Eulerian coordinates 

' pt -b dw{pu) = 0, 

{pu)t + div(p-u 0 n) — pAu — {p + A) Vdivn -b VP(p, 6) = (curlP) x B, 

Cu[{p(^)t + div(pM6')] — kA 9 + 9doP{p,9)divu = 2p\D{u)\‘^ + A(divM)^ -b z/jcurlP]^, ^ ^ 

Bt — curl(M X B) = z/AP, divP = 0, 

where {x, t) G xM’*'. Here the unknown functions p,u = {ui, U 2 , u^Y^, 9 and B = (Pi, P 2 , B^Y^ 
represent the fluid density, velocity, absolute temperature, and magnetic field respectively; D{u) 
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is the deformation tensor and defined by 

D(u) = i [V« + (V«)‘'] . 

The pressure function P{p, 9) is smooth and satishes -Pp(l, 1 ) > 0 and Pe{^i 1 ) > 0 in a neighbor¬ 
hood of (1,1). The constants p and A are the hrst and second viscosity coefficients respectively 
and satisfy the physical restrictions 


p>0, 2/i + 3A>0. 

Positive constants c„, k,, and z/ are respectively the heat capacity, the ratio of the heat conductivity 
coefficient over the heat capacity, and the magnetic diffusivity acting as a magnetic diffusion 
coefficient of the magnetic held. For the sake of simplicity, we assume Cu,Pp{l,l) and Pe(l, 1) to 
be 1. To complete the system fll.ip . the initial data are given by 

(p,M,6»,P)(x,f)|j^o = (Po(a;),Mo(a;),6'o(a;),Po(a;)). (1.2) 

Furthermore, as the spatial variable tends to inhnity, we assume 

lim (po - 1, Mo, 6*0 - l,.Bo)(a;) = 0. (1.3) 

|x|—>-oo 

The compressible MHD systems are combination of the compressible Navier- Stokes equations 
of huid dynamics and Maxwell’s equations of electromagnetism. On the other hand, although 
the electric held E does not apper in fll.ll) . it can be written in terms of the magnetic held and 
the velocity as follows 

E = z/curlP — u X B 

by the moving conductive how in the magnetic held. Obviously, the compressible MHD sys¬ 
tems reduce to the full compressible Navier-Stokes equations when there is no electro-magnetic 
ehect(i.e. B = 0). 

In this paper, we are concerned with the optimal decay rates for higher order spatial deriva¬ 
tives of solutions to the full compressible MHD equations in three-dimensional whole space. Since 
the study of the asymptotic behavior of the MHD equations kept in step with the Navier-Stokes 
equations, we recall some studies on the convergence rates for the compressible Navier-Stokes 
equations with or without external forces. When there is no external force, the convergence 
rates of solutions for the compressible Navier-Stokes equations to the steady state have been 
investigated extensively. First, Matsumura and Nishida [l(] established global existence of small 
solutions in P^-norm and proved that the hrst order spatial derivatives of solutions in Ff^-norm 
converges to zero as the time goes to inhnity in three-dimensional whole space. At the same 
time, Matsumura and Nishida obtained the following convergence rate for all t > 0, 

\\{p-l,u,9- l)(t)||rf 2 < (1 + f)-t, 

if the small initial disturbance belongs to iF^(M^) n L^(M^). For the small initial perturbation 
belongs to only, Matsumura took weighted energy method to show the time decay rates 

||V^(p-l,u,0-l)(f)|U2<(l + f)-t 
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for k = 1,2, and 

\\{p-l,u,e-l){t)\\L^ < (1 + f)-i 

For the same system, Ponce gave the optimal U’ convergence rate 

II V'(p -1,U,9- l){t)\\Lr < (1 + 

for 2 < p < cx) and / = 0,1, 2, if the small initial disturbance belongs to n with 

the integer s > [n/2] + 3 and the space dimension n = 2 or 3. In order to establish optimal decay 
rates for higher order spatial derivatives of solutions, Guo and Wang developed a general 
energy method to build the time convergence rates as follows 


for 0 < / < iV — 1 by assuming the initial perturbation are bounded in H~^{s G [0, |))-norm 
instead of L^-norm. On the other hand, the study of large-time behavior in L^(l < p < oo) 
spaces and pointwise estimates were developed in jbl-lsl. For example, Hoff and Zumbrun p 
studied the isentropic viscous fluid in MF{n > 2) and obtained 


\\{p - 1, pu){t)\\LP 


< 


_IL M — i 
t p 


2 < p < oo, 

l<p<2, 


for all large t > 0, if the small initial disturbance belongs to 0 with the integer 

s > [n/2] + 3, where Ln{t) equals logG-|- t) if n = 2 and 1 otherwise. This result was later 
generalized by Kobayashi and Shibata [oj] and Kagei and Kobayashi 10|, [lit] to the viscous and 


-norm 


heat-conductive fluid and also to the half space problem but without the smallness of L} 
of the initial disturbance. When there is an external potential force F = — V<h(x), there are 
also some results on the convergence rate for solutions to the compressible viscous Navier-Stokes 
equations. For this case, when the initial perturbation is not assumed in L}, the analysis only on 


the Sobolev space yields a slower (than optimal) decay 12, ll3|. If the initial perturbation 

belongs to additionally, Duan et al. [l3, [l^ established optimal decay rates for the solutions 
and its hrst order spatial derivatives as follows 


\w\p-i,u,e 




where k = 0,1. 

Motivated by the study of optimal decay rates for Navier-Stokes equations, the investigation 
of time convergence rates of solutions to the MHD equations has aroused many researchers’ 
interests. First of all, under the il^-framework, Li and Yu 16| and Chen and Tan [l^ not only 
established the global existence of classical solutions, but also obtained the time decay rates for 
the three-dimensional compressible MHD equations by assuming the initial data belong to 
and L‘^{q G [l, |)) respectively. More precisely, Chen and Tan [l7| built the time decay rates 


|V^(p-l,n,H)(t)||H3-. <(l + t) 


. 3 ( 1 . 

2\q 


(1.4) 
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where k = 0,1. These decay rates fll.4l) have also been established by Li and Yu [1^ for the case 
q = 1. Motivated by the work of Guo and Wang j^, Tan and Wang [isl established the optimal 
time decay rates for the higher order spatial derivatives of solutions if the initial perturbation 
belongs to fl H~^ (iV > 3, s G [O, |)). More precisely, they built the following time decay 
rates 




k+s 


where k = 0,1,N — 1. Following the spirit of work [19|, we (see [201 
time decay rates for all t > T*{T* is a positive constant) 




3+2m 


establish the following 


(1,5) 


where m = 2,3. It is easy to see that the time decay rates (II.hh is better than decay rates (11.41) 
since fll.Sp provides faster time decay rates for the higher order spatial derivatives of solutions. 
For the full compressible MHD equations fll.ip . Pu and Guo 2]| established time decay rates for 


the classical solutions in three-dimensional whole space as follows 


- l,u,9 - l,B){t)\\HS-^ < {1 + t) 


.3fl. 

2\q 


( 1 , 6 ) 


where k = 0,1, if the initial perturbation are small in iL^-norm and bounded in [q G [l, |))- 
norm. 

In this paper, we are concerned with optimal decay rates for the higher order spatial deriva¬ 
tives of classical solutions in L^-norm to the full compressible MHD equations in three dimensional 
whole space. For the classical incompressible Navier-Stokes equations, Schonbek and Wiegner 


22 I. |23| applied the inductive argument and Fourier splitting method (see [2^) to establish 


optimal decay rates for higher order derivatives norm after having the optimal decay rates of 
solutions and its first order spatial derivatives at hand. Motivated by 22|, l23| , we move the non¬ 
linear terms to the right hand side of (ll.ip /i and deal with the nonlinear terms as external force 
with the property on fast time decay rates. Then, the application of Fourier splitting method 
helps us to establish optimal time decay rates for the higher order spatial derivatives of mag¬ 
netic field in L^-norm. Since the equation fll.ip is hyperbolic-parabolic type, then the optimal 
decay rate for the second order spatial derivatives of solutions are somewhat complicated. More 
precisely, denoting g = p — 1 and a = 6 — 1, the fll.ip transforms into the system fl2.ip . Then, 
for the homogeneous system fl2.ip fi.e. Si = S2 = S^, = = 0), it is easy to establish following 

energy inequality 

4||v‘(g,«,a)|ii. + c||v'=+‘(t<,ff)||'i. < 0 , (1.7) 

where k = 2,3. In order to apply the Fourier splitting method to build optimal decay rate for 
the second order derivatives norm of solution, we need to rediscover the dissipative estimates for 
g. Then, it is easy to verify the following differential inequality 

^ ■ V^gdx + C\\V^g\\l2 < + || WHi^. 
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The combination of fll.7p and fll.Sp yields directly 

jfi{t) + c[\\v^e\\h + llv“(«,<T)f„0 < 0, (1.9) 

where is equivalent to || V^(^, u, u)||^i. Following the spirit of the Fourier splitting method, 
we need to transform the inequality fll.91) as the form 

jSUt) + ^(2||V“£,|li. + l|V=‘(«,ff)|lP) < 0. 

Then, it is easy to establish optimal decay rates for the second order spatial derivatives of 
solutions in L^-norm(see Lemma (2.61 for detail). Finally, one also establishes time convergence 
rates for the mixed space-time derivatives of solutions. 

Notation: In this paper, we use G M) to denote the usual Sobolev spaces with 

norm || ■ and < p < oo) to denote the usual spaces with norm || ■ \\lv. We dehne 

= {d^Vi\ |a| = /c, i = 1,2,3}, v = (ui,U 2 ,'y 3 )- 

One also denotes the Fourier transform by <^{f) '■= f ■ The notation a < b means that a < Cb 
for a universal constant O > 0 independent of time t. The notation a k, b means a < b and 
b < a. For the sake of simplicity, we write J fdx := fdx. 

The main results on optimal decay rates for higher order spatial derivatives of solutions to 
the full compressible MHD equations fll.ll) - fll.3p can be stated as follows. 

Theorem 1.1. Assume that the initial data (po — l,uo,do — ^ Fl L'^(q G [l, |)) and 

there exists a small constant 5o > 0 such that 


II (PO ~ I 5 '*^05 Sq — 1, Bq)\\jj3 < (5o, 

then the global classical solution (p, u,9,B) of Cauchy problem fll.ip - 
rates for all t > T*{T* is a positive constant), 

\\v\p i)(t)||Hi < c{i + 


-1 




( 1 , 10 ) 

has the following decay 

( 1 . 11 ) 


where k = 2,3. 


Remark 1.1. Obviously, (II.lip provides faster time decay rates for the higher order spatial 
derivatives of classical solutions than fll.bp . Hence, the results in Theorem li.il improve the work 
of Pu and Guo 'Ml- 

Remark 1.2. By virtue of the Sobolev inequality and the results in Theorem li.iL the global 
solution {p,u,9,B) of problem fll.ip - fll.3p has the time decay rates 

||(p- l,u,9 -l){t)\\LP < (^(1+ 

\\V^B{t)\\LP < C'(l + i)-K|-i)-|, 


where 2 < p < 00 , and k = 0,1. 
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Remark 1.3. Although we only established the time decay rates under the H^-framework in 
Theorem \1.1[ the method here can he applied to the H^{N > 3)-framework just following the 
idea as Gao et al. fl£j. Hence, if {po — 1, Uq, Oq — 1, Bq) G n L‘^{N > 3, g G [l, |)) and there 
exists a small constant > 0 such that 


II (Po “ 1) Uq, 6*0 — 1, < ^0) 


then the global classical solutions have the time decay rates 

II -l,u,9- l)(6)||^,iv-. < C{1 + 
WV'^Bit^H^-m < ^(1 

where fc = 0,1,— 1, and m = 0,1, 2,..., N. 


Finally, we establish decay rates for the mixed space-time derivatives of solutions to the 
Cauchy problem fll.ip - fll.3p . 


Theorem 1.2. Under all the assumptions in Theorem li.il then the global classical solution 
{p,u,9,B) of Cauchy problem fll.ip - fll.3p has the time decay rates 


3/1 1 N h I 1 

VVtWII//^-^ + \\^’^Mt)\\L2 + \\V%{t)\\L2 < C{1 + t)-A-.-A-—, 

, 3 / 1 1 ^ k-\-2 


( 1 , 12 ) 


where k = 0,1. 


Remark 1.4. By virtue of the Sobolev ineguality and the results in Theorem li.il the global 
solution {p,u,9,B) of problem fll.ip - (ll.3p has the time decay rates 

\\{pt,Ut,9t){t)\\LP < C(1 
\\Bt(t)\\]^p < C{11) 2(5 p) 


where 2 < p < 6. 

The rest of this paper is organized as follows. In section HI we establish the optimal time 
decay rates for the higher order spatial derivatives of global classical solutions. In section |3l one 
hopes to build the time convergence rates for the mixed space-time derivatives of solutions. 


2 Proof of Theorem 11.11 

In this section, we will establish the optimal time decay rates for the higher order spatial 
derivatives of solutions. By computing directly, it is easy to deduce 

(curlR) xB = {B- V)B - 

and 

curl(M X B) = u{divB) — {u ■ V)B {B ■ V)m — B{divu). 
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Denoting g = p — 1 and a = 6 — 1, the original full MHD equations fll.ip can be rewritten in the 
perturbation form as follows 


Qt + divM = Si, 

ut — pAu — (/i + A)VdivM + Vp + Va = S 2 , 
at — kAct + divM = S 3 , 

Bt — vAB = S' 4 , divH = 0. 

Here Si{i = 1, 2, 3,4) are dehned as 
' Si = — pdivM — u - V Q, 

S 2 = — u ■ Vu — h{Q)[pAu + (/i + A)VdivM] — E{q, a)V q — F{q, cr)Vcr 


+ 9{q) 


B -VB - -V{\B\^) 


( 2 . 1 ) 


( 2 . 2 ) 


(2.3) 


S 3 = — u ■ Vcr — h{Q)(K,Aa) — G{g, cr)divM + g(g) [2p\D{u)\‘^ + A(divM)^] , 

+ g{g){iy\cuT\B\‘^), 

^ S '4 = — M • VH + B ■ Vu — BGrvu, 
where the functions of g and a are dehned as 

P+1 P+1 p+1 

I Pg{0+^,O'+1) T?( ^ Pu{Q+^,(y+P) ^ 

P(P,(^) = -^-1, PiP,<^) = -^-1- 

P + 1 P+1 

To complete the system fl2.ip . the initial data are given by 

(p,M,(T, H)(a;,f)|j^o = (po,Mo,o-o,5o)(a;) (0,0,0,0) as |a;|-)■ cx). (2.4) 

First of all, Pu and Guo (see on Page 521 ) have established the following estimates 

||(p,M,Cr,H)||^^3 < G||(po,Mo,Cro,Ho)||H3, (2.5) 

which, together with the smallness of (5o(see (11.1011 ) and Sobolev inequality, yields directly 

2 ~ ~ 2 

Hence, we immediately have 

|p(p)|<C', |h(p)| < GIpI and |G(p, a)|, |G(p, a)|, |F(p, cr)| < G(|p| + |cr|), (2.6) 


which will be used frequently to derive the temporal decay rates. Furthermore, it is also easy to 
deduce 

|VG(p,a)| = |Gg(p,a)Vp + G^(p,a)Va| < |Vp| + |Va|, 

|V^G(p,a)| = |V(Gg(p,a)Vp +G^(p,cr)Vcr)| 

= |(G^j,(p, cr)Vp + a)Va)Vg + B^(g, a)V'^g 

+ (Eag(g, o-}Vg + E^a(g, a}Va)Va + E^(g, a) V^cr| 

<|Vpp + |Vap + |V2p| + |VV|. 
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In the same manner, we also obtain 


|V%)|,|V^7(^)|<|V^|, 

\^^h{g)\,\V^g{g)\<\Vg\^ + \V^g\, 

\VE{g,a)\,\VF{g,a)\, |VG(^,a)| < |V^| + |Va|, 

\V^E{g, cr)|, \V^E{g, cr)|, |V2G(^, a)| < + |Vcrp + + |VV|. 


(2.7) 


25 


We state the classical Sobolev interpolation of the Gagliardo-Nirenberg inequality, refer to 


Lemma 2.1. Let 0 < m,a < I and the function f G then we have 

l|V“/|U. < ||V™/||i7^||V7irL., (2.8) 

where 0 < 6* < 1 and a satisfy 

First of all, we establish the optimal time decay rates for the second order spatial derivatives 
of magnetic held. 


Lemma 2.2. Under the assumptions of Theorem M . li the magnetic field has the time decay rates 
for all t >0, 

\\V‘^B{t)\\l^ < 0(1 + f)-=^(i-^)-2. (2.9) 


Proof. Taking second order spatial derivatives to fl2.ip multiplying the resultant identity by 
and integrating over then we have 

1 d 


( 2 . 10 ) 


-21 / \V^B\^dx + u / \\/^B\^dx 

2dtj ' ' y ' ' 

= - J V\u-\/B)\/^Bdx + j V^B ■\/u)V^Bdx 

- j V^{BdiYu)V^Bdx 

= h + h + h- 

The application of decay rate fll.Op . Holder, Sobolev and Cauchy inequalities yields immediately 

U = J (VuVH + uV^B)V^Bdx 

< (liVnIUsIlVHlUe + \\u\\l4^^B\\ls)\\V^B\\l2 


< 


\\Vu\\U\^^B\\l,+e\\V^B\\ 


( 2 . 11 ) 


L2 


< (1 + f)"t+^(l + f)"f+^ + e\\V^B\\l2 


< 

r\j 


(1 + f)-^ +£||V'H|| 


L2- 
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In the same manner, it is easy to dednce 


h = - j (V 5 Vm + 

< (||Vn|U3||Vi?|U6 + \\B\\l4V\\\l^)\\V^B\\l. 

< \\\/u\\U\^^B\\l, + \\\/B\\l4VMm+4^^B\\l 
<\\Vu\\U\^B\\l,+e\\V^B\\l, 

<{l + t)-4^i + t)-l4+e\\V^B\\l, 


and 

/3<IIv«||?,.||vb|| 1,.+E||v’B||i. 
<(l+«)-^+£||V’B||i=, 

Substituting fl2.1ip - (l2.13j) into fl2.10p and choosing e small enough, we get 


d 

dt 


V^Bl^dx + V 


V^Bfdx< (1 + t)"^. 


( 2 . 12 ) 


(2.13) 


(2.14) 


Taking third order spatial derivatives to (12.ip ^. multiplying the resulting identity by and 
integrating over then we have 


2 dt 


^'^BVdxFv I W^Bl'^dx 


V^(m-V 5)V^5dx + J V^{B -V^V^Bdx 

{B<\iwu)V‘^ Bdx 


= Ih + Ih + Ih- 


(2.15) 


By virtue of 02.51) . Holder, Sobolev and Cauchy inequalities, it arrives at 

III = J {V4VB + 2VuV^B + uV^B)V^Bdx 

< (llV^nlUallVHlUe + llVnlUallV^HlUe + \\u\\L3\\V^B\\Le)\\V^B\\L2 

< llV^nll^.llV^HlIi. + llVnll^.llV^HlIi. + (e + d)\\V^B\\l,. 
Similarly, it is easy to deduce 

1/2 = - y (V^HVn + 2VBV\ + BV^u)V^Bdx 

< {\\V^B\\L4Vu\\L^ + \\VB\\l4VMl^ + \\B\\l4\^Ml4\\^"B\\l2 

< (llV^'HlU^llVnlUs + \\V^B\\l2\\V\\\l^ + \\VB\\h4^4\\l24V^B\\l2 

< \\V^B\\l4^urH, + ||Vi?||^,||V=^n||i. +£||V^i?||i„ 

and 

Ih < l|V^B||'|„||V«|||„ + IIVBII'I,, II V=>«|||, +£11^5111,. 


(2,16) 


(2.17) 


(2,18) 
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A 

dt 


Substituting fl2.16p - fl2.18p into fl2.15p . then we obtain 

|V=>Bpdx + >. j + (2.19) 

which, together with the decay rates fll.6l) . yields directly 

\V^B\‘^dx + iy J\V^B\‘^dx <{l + ty^. ( 2 . 20 ) 

Adding fl2.15p to fl2.20p . it arrives at 

^ f (|V^5|2 + iV^B^dx + 1^1 (|V^5|2 + \V^Bydx < (1 + (2.21) 


d 

dt 


For some constant R dehned below, denoting the time sphere Sq (see Schonbek [2^) by 


^0 := U e 


then we have 


\v^B\^dx> / 

Jr^/So 


> 


> 


1 +1 
R 

1 +1 


( ^ I tl2| d|2. 


\V^B\^dx 


\ 1 +1 
/ R 



mB\^d^ 


So 




1 + t 



\VB\‘^dx, 


or equivalently 


\V^B\^dx > [\V‘^B\‘^dx - \VB\‘^dx. 


1 + t 


\ 1 + t 



( 2 . 22 ) 


Similarly, it is easy to deduce 

j\V^B\'^dx> j\V^B\'^dx- / \^^B\'^dx. (2.23) 

Substituting fl2.22p and fl2.23p into fl2.21l) and applying the time decay rates fll.6l) . we have 

4 [i\\/^B\‘^+ \V^Bydx+[{\\/^B\^ + iW^B^dx 

dt J l + t J 

< j (ivsr + |v^ii|^)dx + (1 +«)-(;-') 

< (1 + t)-\l + t)"t + 5 + (1 + t)-(l-l) 

< (l+t)"(t + 5) +(l+t)-(f-l). 
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By virtue of g G [l, |), then it is easy to see that 


5-1 

q 


Thus, we have the following estimates 


q 2j 2q - 


d 


Rv 




If choosing 


in fl2.25p . then we get 


R = 


3 + q 
qv 


3+g 


A 

dt 


Multiplying fl2.27p by (1 +1) , it arrives at 

which, integrating over [0,t], gives 

\\V^Bit)\\lr < (1 + [||V'Eoll^i + ^(1 + t)^ 

or equivalently 

Therefore, we complete the proof of the lemma. 


(2.24) 


+ \V^B\‘^)dx < (1 + f)"(l+t). (2.25) 


^ f {\V^B\^ + \V^B\^)dx + f [{\V^B\^ + \V^B\^)dx < (1 + t)-(f+i). 

dt J q{l + t)J 


(2.26) 


(2.27) 


□ 


Next, we establish the following differential inequality for the second order spatial derivatives 
of solutions. 


Lemma 2.3. Under all the assumptions in Theorem M.R then we have 
d 


\W\g, u, a)\\l + (/ill + «:|| VV||i.) < £|| V=^p||i. + (1 + t)"^. (2.28) 


dt 


Proof. Taking k—th.{k = 2,3) spatial derivatives on both hand sides of (I2.ip i . (l2.2p o and (I2.2p i. 
multiplying the resultant identity by V^g, and respectively and integrating over M^, 
then we have 


1 d 

2 dt 


(|V^p|2 + |V^n|2 + iVVHda; + / (/i| V''+'n|" + (/i + A)|V''divn|" + K\V^^^a\'^)dx 




7k+l\2\ 


(2.29) 


= / ■ V^gdx + / V^S2 • V^udx + / • Wdx. 
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Taking A; = 2 in fl2.29p . it arrives at 

lit + + 


j (/i|V^M|2 + (/i + A)|VMivM|2 + K\V^a\^)dx 


= j V^(—^?divM — u ■ Vg)'V‘^gdx + J V^(—m ■ Vm — h{g)[fj,Au + (/i + A)VdivM])V^M(ia: 

+ J V\-E{g, (T)Vg - F(g, a)Va + g(g)lB ■ VB - ^\/i\B\‘^)])V\dx (2.30) 

+ J V^(—M • VcT — h{g){KAa) — G{g^ (T)divM + g{g) \2ii\D{uI + A(divM)^] )V‘^adx 
+ j V‘^{g{g){i'\cm\B\^))V‘^(Tdx. 

Integrating by part and applying fll.Op . Holder, Sobolev and Cauchy inequalities, we obtain 


J V^(—pdivM — M ■ 

= — j (pV^-u + 2V gVu + uV‘^g)V^gdx 

< + llVpIUellV^lUa + ||«|Ua||V2p|U3)||V=^p|U. (2.31) 

< + l|V^p||i.||Vn||^i + \\Vu\\l4V^g\\l,+e\\V^g\\l. 

<(l + t)-M(l + t)-f+U£||V=^p||i. 

Following the idea as fl2.3ip . it is easy to deduce 

J V^(-u ■ Vu)V\dx < II V«|||,.|| V^kII'I,, + £|| V\||i, < (1 + i)-"? + ell (2.32) 

Integrating by part and applying fl2.5p - fl2.7l) . Holder and Sobolev inequalities, it arrives at 


j V^(—h(p)[/iAM + (/i + A)VdivM])V^M(ia: 

< (l|Vfc(e)|U.||V"u||i. + ||ft(t.)||i»||V=>u||i=) ||V=>'t.|U. P-33) 

<3ol|V’«||i.. 

It follows from the integration by part, fll.Op . fl2.6p . fl2.7p . Holder and Sobolev inequalities that 
j V^[-E{g,a)Vg]V\dx 

< \\VE{g,a)hs\\Vg\\L4^ML^ + \\Eig,a)\\L4\^^g\\L4^ML^ 

< l|V(p,a)||i3||V^p||i2 + \\{g,a)\\U\V^g\\h+s\\V\\\h 

< l|V(p,a)||i3||V2p||i. + \Mg,a)\\U\^^g\\l,+e\\VM\l^ (2.34) 

<\\^iAl\\U^^4h + e\\VM\h 

< 4 +tr'^+e\\VM\h. 
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Similarly, it is easy to deduce 

j V^[-F{g,a)Va]V\dx<\\V{g,a)\\l4V^a\\l,+e\\V\\\l, < {1 + t)-"-^ +e\\V\\\l.. (2.35) 
By integration by part, fll.bp . fl2.6p . fl2.7p . Holder, Sobolev and Cauchy inequalities, we get 

J V^{g{Q){B ■VB))V^udx 

< {\\Vg{g)h4B\\L4^B\\pe + \\g{0)U4\VB\4s\\VB\44\\VML^ 

+ \\g{g)\\L4\B\\L4'^^B\\L3\\VML^ 

< \\V^g\\l4VB\\l4V^B\\l, + llVHlIiallV^HlIi. .2.36) 

+ \\VB\\l4V^B\\ls + e\\V\\\l2 

<\\VB\\l4V^B\\lr+e\\V\\\l^ 

< (1 + t)-|+^(l + t)-f+^ + e\\VM\h 

<a+tr-^+4^Mh- 


In the same manner, we obtain 

I j V4-g{g)V{\VB\^))V\dx < (1 + t)" V + e\\VM\h. (2-37) 

Following the idea as fl2.32p and fl2.33p respectively, it is easy to deduce 

j V4-u ■ Va)V^adx < || Vn||^i || W||^i + £|| WHi^ < (1 + ty^T + £||VV||i 2 , (2.38) 


and 




V^(-h(p)VV)VVda: < (5o||VV 


2 

L2- 


(2.39) 


The application of integration by part, fll.6p . fl2.6p . fl2.7p . Holder, Sobolev and Cauchy inequalities 
yields 


V(G(p, (T)divM)V^cr dx 


<(l|VG(g,ff)||i.||V't.|U. + ||Gteff)||i»||V\||i,)||VV|U= 

:£ l|v(i?.ff)|li3||Wti|li. + ||v(t.,<T)||'J,,||v"u||i3 +£||vV|||, 
<l|v(e,ff)|| 5 ,.||W«lli 3 + £||v=><T||i 3 

< (1 + t) 9 + £||V^cr|||2. 

Integrating by part and applying fll.bp . fl2.6p . fl2.7p . Holder and Sobolev inequalities, we get 


J V^(5'(p) [2fi\D{u)\'^ + A(divM)^] )'V'^adx 
- j v(^(p)|VM|^)vVdx 

< 4^g{g)\\L4^u\\y + ||^?(p)|Uoo||Vn|U3||V2n|U6)||VV|U2 (2.41) 

< \\\/4\\l4^4\\h + l|Vn||^i||V=^n||i2 + £||VV||i2 
<l|V(p,n)||^i||V\||^.+£||VV||i2 

< (1 + t) 1 + £||V^cr||^2. 
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In the same manner, it is easy to dednce 


V‘^{g{Q){u\cnT\B\‘‘^))V^adx < (1 + t) «’ + ^11 V^(j||^2 


(2.42) 


Substituting fl2.3ip - fl2.42p into fl2.30p and applying the smallness of Sq and e, it is easy to deduce 

^ f (IV'pI' + IV'nl' + IVVHdx + f (p|V\|2 + < (1 + t)-V + e\\V^g\\l2. 


Therefore, we complete the proof of the lemma. 


□ 


Furthermore, we establish the following differential inequality for the third order spatial 
derivatives of solutions. 


Lemma 2.4. Under all the assumptions in Theorem \l.l[ then we have 
d 


dt 


\V\q, u, a)||i2 + (/ill + «|| VVlIi^) < (1 + t)-^ + e|| V^plli^. 


(2.43) 


Proof. Taking /c = 3 in fl2.29p . it is easy to see that 


1 d 


J + I+ j + (/i + A)|VMivM|2 + i/|vvndx 

= y V^(—pdivM — u ■ '\/g)'V^Qdx -|- ■ Vu — h{g)[p,Au + {p, + A)VdivM])V^Mda: 

+ j V%-E{g, a)Vg-F{g,a)Va + g{g)[B-VB-^V{\B\^)])V^udx 
+ J V^(—M ■ VcT — h{g){KAa) — G{g, a)divu + g{g) \2pL\D{u)\^ + A(divM)^] )'V^adx 

+ J V^{g{g){h'\cuT\B\‘^))'V^adx. 

By virtue of decay rates fll.bp . fl2.5p . Holder, Sobolev and Cauchy inequalities, we obtain 


(2.44) 


J '^^{—gdivu)V^gdx 

< bllL-||V^ii|U2||V='p|U2 + ||Vp|U3||V3n|Ua||V3p|U2 

+ IIV^pIIlsIIV^m||l6||V^p||l2 -|- ||Vm||loo||V^p||^2 ^-2 

< IIV^pII^iIIV^ ii||i2 + (e + do)(llV=^p||i2 + IIV^ii||i2) 

< (1 + /)-f+^(l + /)-|+^ + (e + do)(||V=^p||i 2 + IIV^ii||i 2 ) 

<(l+t)-'?+(e + do)(||V=^p||i2 + ||V^ii||i2). 

Similarly, it is easy to deduce 

j V^{-u-Vg)V^gdx<{l + t)-"-i^+ {e + do)\\V^g\\h. (2.46) 
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With the help of decay rates fll.6p . fl2.5p . Holder, Sobolev and Cauchy inequalities, we get 


J V^(—M • Vn)V^M dx 

< (||Vm||l3||V\||l6 + \\u\\l3\\'^^u\\l6)\\V%\\l2 

< IIVnll^ill+ (£ + 5o)||V^n||i. (2-47) 

< (1 + + (£ + <5o)||V^nlli^ 

<(l + t)-'^ + (^ + <5o)||V^n||i.. 

By virtue of the integration by part and applying decay rates fll.bp . fl2.5l) - fl2.7p . Holder and 
Sobolev inequalities, we deduce 

J V^(—h(p)[/xAM + (/i + A)VdivM])V^M(ia: 

~ J dx 

< (llVpIliellV^IUe + ||V2p|U3||V2n|U6)||V"n|U2 (2.48) 

+ (||Vp|U3||V3n|U6 + ||h(p)|Uoo||V"n|U2)||V"n|U2 
<||V2p|||.||V=^n||i. + (5o + £)||V^n||i. 

<(l + t) 'J + (e + (5o)|| V"^m||^2. 

The application of decay rates (11.61) . (12.6p . (12.7p . Holder and Sobolev inequalities yields directly 


- J V^[-E{g,a)Vg]V^udx 

= J['V'^E{g, a)Vg+ 2VE{g, cr)V^p + E{g, a)V^g]V'^udx 

< (||Vp + Va||ie||Vp|U6 + ||Vp|U3||V2p + VV|U6)||VS|U2 

+ (||Vp + Va|U3||V2p|U6 + ||E(p,a)|Uoc||V3p|U2)||V"n|U2 (2.49) 

< l|V^(p,cr)||^2 + ||V(p,a)||i3||V^(p,a)||i2 +e||V^n||i2 

+ l|V(p,a)||^,||V=^p||i2 

< l|V^(p,a)||i2 + ||V(p,a)||^.||V3(p,a)||i2 +£||V^n||i2 

< (1 + t) 1 + e:|| V"^'u|||2. 

In the same manner, it is easy to deduce 

j V^[-E{g, a)Va]V^udx < (1 + t)"^ + £|| V^M||i 2 . (2.50) 
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Integrating by part and applying fll.bp . fl2.6p . fl2.7p . Holder and Sobolev inequalities, we get 

J \/^gig){B ■VB))V\dx 

= - j{V‘^g{Q)BVB + 2Vg{Q)V{BVB) + g{Q)V^{BVB))V%dx 


< (IIIV^pI + |VpnUa||H|U6||Vi?|Ua + \\Vgh4^B\\l,)\\V\h. 

+ (llVpIUallHlUellV^HlUe + || || Vi^lUa || || 

+ ||^?(p)|Uoc||i?|U.||V3i?|U2||V^n|U2 

< llVHlIi.llV^HlIi. + llV^HlIi.llV^HlI^, +e||V^n|li. 
<\\VB\\U\^^B\\l,+e\\V\\\l. 

<{i+tr-^+4^M\h- 

Similarly, it is easy to deduce 

i j V’(-9(e) V(| VBp))V’«di <(i + t)-‘-?+ £|| v-Mlli,. 
Following the idea as fl2.47p and fl2.48p respectively, it is easy to deduce 

j V^(-u ■ Va)V'‘adx < ||V(u,<T)||y ||V“(«,<T)||i, + (£ + *)l|VV||i, 

<(l + i)-‘? + (9 + «l|VV||i., 

and 


(2.51) 


(2.52) 


(2.53) 


K j V^(-h(p)VV)VV dx 

< ||V'^9|||„||VV||i. + (9 + i)||VV||i. (S'S-i) 

<(l+()-"? + (^ + ^)l|VV|||.. 

Integrating by part and applying decay rates fll.bp . fl2.6p . fl2.7p . Holder and Sobolev inequalities. 


it arrives at 


J V^[—G(p, (T)divM]V^(J(ia; 

= J [V^G(p, cr)divM + 2VG(p, (T)VdivM + G{q, a) V^divM]V"^crcia: 

< [||Vp + Va||ie||VM|U6 + ||V2p + VV|U3||VM|U6 

+ ||Vp + Va|U3||V\|U6 + ||G(p,a)|Uoo||V3w|U2]||VV|U2 

< l|V^(p,n,a)||i. + ||V(p,a)||^.||V^n||^. +£||VV||i. 

+ l|V(p,a)||^.||V\||i. 

< \\^\^,n,a)\\% + ||V(p,a)||^.||V^n||^. +£||VV||i. 

<(l + t)“^+^l|VV||i.. 


(2.55) 
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Integrating by part and applying decay rates fll.bp . fl2.6l) . fl2.7p . Holder and Sobolev inequalities, 
we get 

J [2r\D{u)\‘^ + \{divuY])'V^adx 

~ [ V\g{g)\Vu\^)V^adx 


= j {V^g{g)\Vu\^ + AVg{g)VuV^u)V^adx 

+ J (2g{g)V^uV‘^u + 2g{g)VuV^u)V^(jdx 

< (||Vp|U^||Vp|U6||Vn||ie + ||V^p|U6||Vn||ie)||VV|U2 

+ (llVpIUellVwlUellV^MlUa + || V'ulUs || V^lUa) || VV|U2 
+ II Vti II L°° II V^mII II V^(t|| 

<l|V^p||^i||V\„,. , . 

<l|V^(p,«)||^i||v2,dl?.+HIVVII^ 

< (1 +t)“V +e||vV| 

In the same manner, it is easy to deduce 


(2.56) 


^ m ||^ i || V ^ m ||^2 +£|| V ^ cr ||^2 


74_||2 

Il2- 


V“(g(t.)(i^|curlBp))VV* < ||V"(g,B)||y||V 2 B|||,. +£||VV||i. 


(2.57) 


< (1 + t) 9 + £||V cr||^2. 


Substituting (I2.45I) - (I2.57I) into (I2.44p and applying the smallness of 5o and £, it is easy to deduce 

^ j (I V"p|2 + I V^'nl' + I VVHdx + j + KlVV^dx < (1 + + £|| V^p||i2. 

Therefore, we complete the proof of the lemma. □ 

Furthermore,, we establish the dissipative estimates for the density. 

Lemma 2.5. Under all the assumptions in Theorem \l.l[ then we have 

^ y ■ V^gdx + j iV^pl^dx < Ci(|| VV||i2 + || V^m||^i) + Ci( 1 + t)"^. (2.58) 

Proof. Taking V^—th spatial derivatives on both hand sidse of fl2.iP o. multiplying by and 
integrating over then we get 

J (VV • + I= j [/iAV\ + (/i + A)VMivM - VV + V^pda:. (2.59) 

In order to deal with the term f V^Ut ■ V^gdx, we turn to time derivatives of velocity to the 
density and apply the transport equation fl2.ip ^. More precisely, we have 

/ V'^Ut ■ V^gdx = 


d 

dt J 

f V^U' 

■ - J 

[ ■ V^gtdx 


d 

dt J 

f 

■ V^gdx + J 

[ V^divM • V'^gtdx 

(2.60) 

d 

dt J 

[ 

■ V^gdx — J 

f V^divM ■ V^(divM + pdivu + uVg)dx 
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Substituting fl2.60p into fl2.59p . it is easy to deduce 

^ y + j\V^Q\^dx 

= J I V^divul^dx + J V^divM • V^(^?divM + 

+ [ • ^^gdx + [ [/iAV^u + (/i + A)VMivu - Va]V^gdx. 


(2.61) 


Integrating by part and applying decay rates fll.6p . Holder and Sobolev inequalities, we obtain 
J V^divn • V^(pdivn + nVp)(ia; 

= — J V^divn • V(pdivn + uVg)dx 

~ I|V'^-u||l2(||VpIIlsIIV- u||i6 + ||p||l6||V^“IU3 + ||m||l6||V^pIIls) 


(2.62) 


< 


l|Vp||^i||V^n||i. + ||Vp||i.||V^n||^. + ||Vn||i.||V^p||^, + £||V^n||i. 


+e\\V\f 


< (1 + ty - -r t|| V u||^ 2 . 

Following the idea as in Lemma 12.41 we deduce immediately 

[ • V^gdx < (1 + t)-V + e\\V^g\\l,. 


(2.63) 


On the other hand, it is easy to see that 


J [pAV^n + (p + A)V^divM - < || VV||i2 + ||V^n||i2 + e\\V^g\\l2. (2.64) 

Plugging (I2.62p - (l2.64p into (I2.6ip . we complete the proof of the lemma. □ 

The optimal decay rates for the second order spatial derivatives of global classical solutions 
are stated in the following lemma. 

Lemma 2.6. Under all the assumptions in Theorem M . li then the global classical solution {g, u, a) 
of Cauchy problem fl2.1l) - fl2.4p has the decay rates 

l|V'p(t)||l,i + \\^Mt)\\li + IIVV(f)||^i < 0(1 + t)-3(i-5)-2 (2.65) 


for all t > T*(T* is a positive constant defined below). 

Proof. Adding (12.281) with (12.431) . it is easy to deduce 

^||V2(p,n,a)||^i + (/^||V='n||i2 + «||VV||^0 < + C2{1 + . (2.66) 

Multiplying (I2.58p by and adding with (I2.66p . we obtain 

j^M^t) + OsdlV='p||i2 + IIV^ll^i + IIVVlI^O < 04(1 + t)-^, (2.67) 
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where 

Mlit) = \\V\g,u,a)\\l, + j V^u-V^gdx. 

Applying the Young inequality and the smallness of Sq, we have the following equivalent relations 

Cg\\V^(Q,u,a)f„, < Mlit) < C5||V=(ti.ti.<T)||y. (2.68) 

From the inequality fl2.67p . it is easy to deduce 

jMKt) + A(|| vvili. + II v=>£.||i, + II + II VV||5,.) 

< jM^t) + A(2||V’i)||i» + 2||V=«|I?,. + 2||VVliy) 

< C*4(l + t) 9 , 


or equivalently, 


+ ^(11 V^p||^2 + IIV^p||^2 + IIV^mII^i + IIV^all^i) < C4(l + t) ^ 


Similar to (12.221) . we have 


\V^u\‘^dx > 


R 


and 


\V%\^dx > 


1 + t 
R 

1 + ^ . 


\V\\^dx-(^—] \Vu\'^dx, 


\V^u\^dx 


[ 

(- 


1 + t 


R 



\ 1 + t 





Adding fl2.70p with fl2.7ip . we obtain 


R 


( R 


iiv'flii,. > “ (ttt) 


In the same manner, it is easy to deduce 


R 


and 


11 ^Vlli. > Yfi\\^‘e\\h 


l|VVf„. > yYiIWII?,. 




\ 1 + t 


1 l|Vp||i2, 




Combining fl2.72l) - fl2.74l) with fl2.69p and applying the decay rates fll.bp . then we get 

R 


d , ,2. . C*3 

W + Y 


1 + t 


(Iiv^p|li2 + ||v^«ll^2 + livvii^0 + liv^p|li2 


(z -(l|Vp||i 2 + ||V-u||^i + ||Vcr||^i) + (1 + t) 


< 

~ \l + t 


< (1 + t)-2(l + t)-U-2) + (1 + 

< (l + t)-(f+5) + (l + t)-(f-l) 


(2.69) 

(2.70) 

(2.71) 

(2.72) 

(2.73) 

(2.74) 


(2.75) 
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where we have used the fact fl2.24p . For some large time t > R — 1, we have 

R 


which implies 


R 


1 +1 




< 1, 


3„l|2 


1 + t 


Plugging (12.7611 into (12.7511 . it is easy to deduce 

C.R 


d , . Cy.sJ 

WYT) 


v^eliy + iiv’ti|iy + IIwiiy) < (i +«) 


which, together with the equivalent relation (I2.68p . gives directly 

Choosing R = in (I2.77p . then we have 




(2.76) 


(2.77) 


(2.78) 


Multiplying (I2.78P by (1 + t) 'j , it arrives at 

<(l + t)-i (2.79) 

The integration of (I2.79p over [0,f] yields 

M2{t)<{l+t) 1 (M|(0) + C'(l + f)2), 

which, together with the equivalent relation (I2.68p . gives 

iivv(i)iiy + iiv=«(i)iiy + iivv(«)f„. < c(i+()-("■"’>, 

for all t > T* := — 1- Therefore, we complete the proof of the lemma. □ 

Finally, we establish optimal decay rates for the third order spatial derivatives of magnetic 
held. 

Lemma 2.7. Under all the assumptions in Theorem M . 1[ then the magnetic field has the following 
time decay rate for all t > T* 


a 

dt 


3+g 


{l+t)^Ml{t) 


V^B{t)\\l2 < C(1 +t)“^(^“^)“l 


(2.80) 
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Proof. Combining the time decay rates fl2.28p . fl2.65p with fl2.38p . we get 

d 


^ / \V^B\‘^dx + u / \V^B\‘^dx 




3,,l|2 


lL2 


Combining fl2.40p . fl2.8ip and the time decay rates fl2.28p . we obtain 

4- [\V^B\^dx + [\V^B\^dx 

dtj ' ' 1 + t / ' ' 


< 


(iT^ 


V^Bfdx + (1 + t) 9 


<(l+t)-U+2j+(l+t)-^ 
By virtue of g G [l, |), it is easy to see that 

t_ C + = s(i^9) 

q \q 2 ) 2q 


> 0 . 


Then, we deduce from the inequality (I2.82p that 

4- [\V^B\^dx + ^^ [\V^B\^dx + 
dtj' ' 1 + tJ' ' ~ ^ ^ 

Choosing R = and multiplying fl2.83p by (1 + 1) 9 , then we have 

|i(i + ()“iiv=>B|iy <(i + ()-i 


Integrating fl2.84p over [0,t], it arrives at 

\\V^B{t)\\l, < (1 + t)-(f+2)(||V=^i?o||i2 + C(1 + f)^), 
which implies the following time decay rate 

\\V^B{t)\\l2 < C(1 +t)"^(9 + 5). 

Therefore, we complete the proof of lemma. 


Proof of Theorem 11.11 With the help of Lemma 12.21 Lemma 12.61 and 
complete the proof of the Theorem 11.11 


(2.81) 


(2.82) 


(2.83) 


(2.84) 


□ 

Lemma 12.71 we 
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3 


Proof of Theorem 


1.2 


In this section, we establish time decay rates for the mixed space-time derivatives of solutions. 

Lemma 3.1. Under the assumptions in Theorem M . 1[ the global solution {g,u,a,B) of problem 
(EIID-dSSD has the time decay rates 

3/1 1 \ fe+1 

+ \\V’^ut{t)\\L^ + ||VVi(t)|U2 < ^(1 , 


where k = 0,1. 


Proof. First of all, applying fl2.ip i . fll.lip . Holder and Sobolev inequalities, it arrives at 

= II - divM - pdivM - M ■ Vp ||^2 

<l|V«||i2 + ||p||ia||V«||i3 + ||«||ie||Vp||i3 

< (1 + t)-l + -^ + (1 + + + ^ 

Similarly, we obtain 

II = II — VdivM — V(pdivM — m ■ Vp )||^2 

<l|V^«||i2 + ||Vp||i3||V«||ia + ||p||i3||V\||ie 

+ l|Vu||ia||Vp||i3 + ||«||ia||V^p||i3 (3.2) 

<l|V^p||^i + ||V^«||^, 

< (1 


v^ftiii2 = 

l|v'(- 

-divM 

— pdiVM — M • Vp)| ^2 


< 

r\j 

I|v^m| 

|2 _L 

Il2 + 

llV^pIliallV^llia + llpIlioollV^ 

n\\h 


+ 1 V^l Lsl 

v^«lli3 + hilioo||v=^p|li2 


< 

r\j 

I|v^m| 

|2 -L 

Il2 + 

I|V^p||^i||V2«||^, + ||V«||^3| 


< 

r\j 

||V\| 

|2 _L 

Il2 + 

IIV^pII^i 



< (1 +f) 5 2 . 

Combining fl3.ip - fl3.3p . then we have the time decay rates 

||V'=Pi(t)||^2-. < ^(1 (3.4) 


where k = 0,1. Secondly, in view of the equation fl2.ip „. fll.lip . fl2.6p . fl2.7p . Holder and Sobolev 
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inequalities, it is easy to deduce 

= ll/iAw + (/i + A)VdivM — Vg — V(j ||^2 
+ II — M ■ Vu — h{g)[fiAu + (/i + A)VdivM] — E{g, ct)V ^||^2 

+ II - Fig, a)Wa + gig)[B ■ VB - ^V(|i?p)]||i 2 

~ + II cr)||^2 + ||m||^6|| VmII^s + ||^(^)|Il°o||V^'u||^2 

+ \mg,a),Fig,a))\\U\Vig,a)\\l, + \\gig)^^^^^^^ 

<\Mn,B)\\l. + \\Vig,a)\\l, 

Similarly, it is easy to deduce 

II VMt||i 2 = ||V(/iAM + (/i + A)VdivM - - Va)||i 2 + ||V52||i2 

< l|V=^«||i. + ||V2(^,a)||i2 + \\Vu\\U^\\\l, 

+ l|V(^,cr)||^i||V2(^,a)||i2 + ||Vi?||^,||V25||^i (3.6) 

<l|V^(«,i?)||^i + ||V2(^,a)||i2 

< il + t)~l~F 

Combining (I3.5p - (l3.6p . then we have the time decay rates 

||V^Mt(t)||i2 < ^(1 + (3.7) 

where fc = 0,1. Furthermore, it is easy to deduce 

ll'^i|lL 2 = ||kA(T — divM — u ■ Va — hig)iKAa) — Gig, cr)divM ||^2 
+ Il5'(^) [2/i|F)(M)p + A(divM)2] + 5((^)(z/|curl5n||^2 
< ||VV||i2 + ||V«||i2 + ll^lliellVallia + ||M^>)llioc||VV||i2 
+ l|G(p,a)||i.||V«||i2 + ||^7(^>)||ioc||V«||i3||Vu||ia (3.8) 

+ ll<?(^)llioc||V5||i3||V5||ie 
<||VV||i2 +II V(u, 5)11^1 

Similarly, it is easy to deduce 

II VcTt||i2 = II V(«:Aa - divM)||i2 + II V^slli^ 

< IIvV||i2 + ||v\|li2 + ||Vk||^i||vV||^3 

+ l|V(^>,a)||^,||V^«||i2 + ||V(p,«)||^3||V^«||^, 

+ l|V(^,i?)||^l||V^i?||^3 ^ ^ 

<\\V\a,u,B)\\l, 

< (i + t)"i"i 
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In view of fl3.8p - fl3.9l) . then we have the time decay rates 

IIVVi(t)||i2 < C{1 + t)-3(|-5)-(^+i)^ 

where /c = 0,1. Finally, it follows from the fl2.iP o. Holder and Sobolev inequalities that 

ll-^illia = ||z^Ai? — u ■ VB + B ■ Vu — i?divM ||^2 

<l|V^H||i. + ||«|lie||Vi?||i3 + ||i?||ia||V«|li3 
<l|V^H||i2 + ||V(n,H)||i.||V(n,5)11^3 

< (1 + t)-hh + (1 + t)-| + i(l + t)-f + 5 

< (1 + 

Similarly, it follows from fl2.iP o. fl2.30p - fl2.32p that 

II Vi?t ||^2 = ||V(i^Ai? — u ■ VB + B ■ Vu — i?divM )||^2 
<||V3H||i2 + ||Vn||^2||Vi?||^2 
<(l + t)-f-f+ (l + t)-f+i 

By view of fl3.1ip and fl3.12p . it is easy to obtain 

\\V^Bt{t)\\l2 < C{1 + 

where /c = 0,1. Therefore, we complete the proof of the lemma. 

Proof of Theorem [list With the help of Lemma [3Tl we complete the proof of the Theorem 

o 
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